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ABSTRACT 

Aims. We aim at investigating with high angular resolution the NH 3 /N 2 H + abundance ratio toward the high-mass star-forming region 
AFGL5142 in order to study whether the NH3/N9IL ratio behaves similarly to the low-mass case, for which the ratio decreases from 
starless cores to cores associated with young stellar objects (YSOs). 

Methods. CARMA was used to observe the 3.2 mm continuum and N 2 H + (1-0) emission toward AFGL5142. We used NH 3 (1,1) and 
(2,2), as well as HCO + (1-0) and H' 3 CO + (1-0) data available from the literature to study the chemical environment. Additionally we 
performed a time-dependent chemical modeling of the region. 

Results. The 3.2 mm continuum emission reveals a dust condensation of ~ 23 Mq associated with the massive YSOs, deeply 
embedded in the strongest NH3 core (hereafter central core). The dense gas emission traced by N 2 H + reveals two main cores, the 
western core of ~ 0.08 pc in size and the eastern core of ~ 0.09 pc, surrounded by a more extended and complex structure of ~ 0.5 pc, 
mimicking the morphology of the NH3 emission. The two cores are located to the west and to the east of the 3.2 mm dust condensation. 
Toward the central core the N 2 H + emission drops significantly, indicating a clear chemical differentiation in the region. The N 2 H + 
column density in the central core is one order of magnitude lower than in the western and eastern cores. Furthermore, we found 
low values of the NH3/N2IL abundance ratio ~50-100 toward the western and eastern cores, and high values up to 1000 associated 
with the central core. The chemical model used to explain the differences seen in the NH 3 /N 2 H + ratio indicates that density, and in 
particular temperature, are key parameters in determining the abundances of both NH3 and N 2 H + . The high density (n = 10 6 cnT 3 ) and 
temperature (T a 70 K) reached in the central core allow molecules such as CO to evaporate from grain mantles. The CO desorption 
causes a significant destruction of N 2 H + , which favors the formation of HCO + . This result is supported by our observations, which 
show that N 2 H + and HCO + are anticorrelated in the central core. The observed values of the NH 3 /N 2 H + ratio in the central core can be 
reproduced by our model for times / = 4.5 - 5.3 x 10 5 yr while in the western and eastern cores the NH3/N 2 H + ratio can be reproduced 
by our model for times in the range 10 4 - 3 x 10 6 yr. 

Conclusions. The NH3/N 2 H + abundance ratio in AFGL5142 does not follow the same trend as in regions of low-mass star formation 
mainly due to the high temperature reached in hot cores. 

Key words, astrochemistry - stars: formation - ISM: individual objects: AFGL5142 -ISM: clouds -ISM: molecules -ISM: abun- 
dances 

1 . Introduction the same trend, with a high NH 3 /N 2 H + ratio, up to 300, in the 

cores with starless properties, and around 50 for the cores as- 
£ ; It is well known that N-bearing molecules, such as NH 3 and socia ted with YSOs, clearly showing that chemical differentia- 
. N 2 H + , are excellent tracers of the interstellar dense gas because tion is i mport ant in the region. All these studies show that the 
j_ ■ none of these mole cules deplete onto dust grains until densities NH 3 /N 2 H + abundance ratio is consistent with being a 'chemical 
03 ' reach ~ 10 6 crn^ feerein & Langet1ll997l: rTafalla et al J 120041: c i oc k'. 

| FloweretaU200 l. Thus observations of these dense gas tracers ^ In order to investigate the behavior of the NH 3 /N 2 H + abun- 
have become a powerful tool to study the sites of star forma- dance ratio - m high . mass star .f orming reg ions, and see if the 
tion. However, only a few observational studies have focused N h 3 /N 2 H + ratio in massive YSOs shows the same trend as in 
on the comparison of NH 3 and N 2 H + cores, and these studies, the low . mass ^ we carried out observations and chemical 
which were carried out toward low-mass star-forming regions, modeling of the high . mass s tar-forming region AFGL 5 142, lo- 
find that the NH 3 /N 2 H + abundance ratio is around 60-90 close cated at a difjtance of Lg kpc (Snell et al 198g) in the Perseus 
to the young stellar objects (YSO^yJiile it rises up to 140- am The region was selected as good candidate t0 study the 
190 in starless cores ( |Caselli et al. | l2002at i HotzeletalJ |2004t N h 3 /N 2 H + ratio because the high angular resolution NH 3 emis- 
i Fnesen et al. || 2010|) . In the intermediate-mass cores surrounding sion reyeals seyeral dense cm which haye very different tem _ 
the high-mass star IRAS 20293+3952, | Palau et alj <|2007b find peratures> and one of them is associate d with star formation in 
clustered mode with the presence of hot cores. The region con- 
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dHunteret all 1 19951) IRAS 05274+3345, with a bolometric lu- 
minosity of 3.8xl0 3 Lq (ICarpenter et alJll990h. li es near the 
western edge of the NH3 gas dEstalella et al.lfl993l) and a lack 
of dense gas emission associated with it suggests that the IRAS 
source is a more evolved region , which is consiste nt with the de- 
tection of optical nebulosity bv lEiroa et al. ( 1994]). Bas ed on its 
position in the J vs. J - H diagram, Chen et al . (2005) classify 
the IRAS source as a B2 star of 3 Mvrs. iTorrelles et alj d 19921) 
detect a faint radio continuum source, IRAS 05274+3345 East 
(hereafter referred as AFGL 5142), about 30" to the east of 
the IRAS source, which coincides with the p eak position of the 
dust emissi on at 3.4 mm ([Hunter et alj |1999) and the NH3 emis- 
sion peak dEstalella et alJH993l) . Assuming that the centimeter 
emission arises from an optically thin H 11 region, the flux den- 
sity is equivalent to a zero-a ge main sequence (ZAMS) star of 
spectr al type B2 or earlier ( Torr e lies et al.l [l992l : iHunter et alJ 
1995b with a lumin osity of 4 X 10 3 Lq. In the near-infrared, 
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Hunt er et al.l ([T995) report a cluster of 28 embedded sources in 
a region of 1' (~0.3 pc in radius ) near the positio n of the ra- 
dio continuum source. Recently, lOiu et all d2008l) carried out 
Spitzer Space Telescope IRAC and MIPS observations toward 
AFGL 5142. The authors identify 44 YSOs, 20 of them clus- 
tered around the central massive star, and surrounded by a 
more extended and sparse distribution of young stars and proto- 
stars. In addition, high angular resolution observations of maser 
emission revea l a cluster of H2O and CH3OH masers in an 
area of ~ 5" dHunter et all fl995t iGoddi & Moscadellil 12005 
iGoddi et all2007 ). Subsequent observations with the Very Large 
Array (VLA) in the A configuration at 8.4 GHz and with the 
Submillimter Array (SMA) a t 1 .2 mm wi t h ~ 1 " of angular res- 
olution were carried out by Zha ng et all (|2007), who find that 
the centimeter source is reso lved into three peak s, and that the 
3.4 mm source detected by IHunter et al.l (l999) actually con- 
sists of five millimeter sources. Regarding the molecular outflow 
emission in this region, Zhan g et al.l d2007l) identify three molec- 
ular outflows in CO (2-1) and SO 65-54, one of them coinciding 
wi th an HCO + outflow and the well-collimated SiO jet detected 
by IHunter et al.l d!999l) . All three outflows appear to originate 
from the dust condensation in a region of about 3". All this in- 
formation indicates that active star formation in clustered mode 
is taking place in AFGL 5 142. 

The dense gas emission in this region has been studied with 
single-di sh telescopes in NH3, CS HCN HCO + , CH3OH, and 
CTLCN dVerdes-Montenegro et alj|1989t lEstalella et alj [1991 
IHunter etal.Hl 9951 119991: ICesaroni et all 19991) . In particular, the 
dense gas emission traced by the NH3 molecule has been ob- 
serve d with high angular resolution using the VLA ( Zhang et alj 
|2002|) . The high angular resolution NH3 emission consists of 
a central and compact core associated with the dust condensa- 
tion, harboring at least three intermediate/high-mass young stars, 
surrounded by fainter NH3 cores located in a more extended 
structure with no signs of stellar activity (no maser nor molecu- 
lar outflow emission associated) indicating that the region may 
harbor cores in different evolutionary stages. The presence of 
cores containing massive star(s) together with cores with no star- 
formation activity makes this region a good choice to study how 
the NH3/N2H + ratio behaves in high-mass star-forming regions. 

In this paper we report CARMA observations of the contin- 
uum emission at 3.2 mm and the dense gas traced by N2H + (1- 
0) toward AFGL 5 142. The paper layout is as follows: in § 2 we 
summarize our observations and the data reduction process. In 
§ 3 we present the main results for the continuum and N2H + 
molecular line emission. In § 4 we analyze the molecular emis- 
sion of several species by computing their column density maps. 
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Fig. 1. 3.2 mm continuum emission toward AFGL 5142. 
Contour levels range from 3 to 60 cr in steps of 3 cr, where 
cr is the rms of the map, 0.7 mJybeairT 1 . The synthesized 
beam, shown in the bottom right corner, is 2'.'53 x 1'.'71, RA.= 
68.8°. White crosse s indicate the millimeter sources detected by 
IZhang et all d2007l) at 1 .2 mm with the SMA. 



In § 5 we show the main results of the chemical model to qual- 
itatively reproduce the abundances of the region. Finally, in § 6 
we discuss our findings, and we list the main conclusions in § 7. 



2. Observations 

The Combined Array for Research in Millimeter-wave 
Astronom >Q (CARMA) was used to observe the 3.2 mm 
continuum and the N2H + (l-0) emission toward AFGL 5142. 
CARMA consists of six 10 m and nine 6 m antennas located at 
2200 meters elevation at Cedar Flat in the Inyo Mountains of 
California. The observations were carried out on 2007 February 
4 and March 11 using the array in the C configuration with 
14 antennas in the array. The projected baselines ranged from 
26 to 370 m. The phase center was set at a = 05 h 30 m 48 s .02, 
6 = +33°47'54'.'47. The FWHM of the primary beam at the 
frequency of the observations was 132" for the 6 m antennas 
and 77" for the 10 m antennas. System temperatures were 
around 250 K during both days. 

The digital correlator was configured to observe simulta- 
neously the continuum emission and the N2H + (l-0) group 
of the hyperfine transitions (93.176331 GHz, in the lower 
sideband). The continuum data were recorded in two ~500 MHz 
bands covering the frequency ranges 93.44-93.88 GHz and 
96.74-97.18 GHz from the receiver lower and upper sidebands, 
respectively. We used two consecutive bands of 8 MHz of 
bandwidth with 63 channels in each band, providing a spectral 



1 Support for CARMA construction was derived from the Gordon 
and Betty Moore Foundation, the Kenneth T. and Eileen L. Norris 
Foundation, the Associates of the California Institute of Technology, the 
states of California, Illinois, and Maryland, and the National Science 
Foundation. Ongoing CARMA development and operations are sup- 
ported by the National Science Foundation under a cooperative agree- 
ment, and by the CARMA partner universities. 
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Table 1. Parameters of the hyperfine fits to the N2H + (1-0) line for the central position of each core 
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Position" 


A< 


V LSR C 


Av 




T 

' ex 


2V(N 2 H + ) 


Core 


a(J2000) 


(5(J2000) 


(Kkms- 1 ) 


(kms -1 ) 


(km s" 1 ) 


r d 

' m 


(K) 


(xlO 12 cm" 2 ) 


central 


05:30:48.02 


33:47:53.4 


0.51 ±0.04 


-0.97±0.06 


1.1 1±0.10 


0.10±0.03 


8.2±1.6 


2.8±0.9 


eastern 


05:30:48.73 


33:47:52.9 


1.35±0.02 


-2.92±0.01 


1.60±0.02 


0.10±0.02 


16.7+1.1 


18.5+1.2 


western 


05:30:47.47 


33:47:51.7 


2.45±0.01 


-3.57±0.01 


1.31+0.03 


0.27±0.02 


12.4±0.8 


14.4+1.2 



Notes. a Positions, given in (h m s) and (° ' "). b A = f{J v {T a ) — J v (T bg )), where / is the filling factor (assumed to be 1), T ex is the 
excitation temperature of the transition, 7b 6 is the background radiation temperature, and J V {T) is the intensity in units of temperature, 
J V (T) = (hv/k)/(exp(hv/kT) - 1). c u LSR corresponding to the F t F = 01 — > 12 hyperfine component . d Optical depth of the main line, 
F\F = 23 — > 12, obtained from the fit. e Position of the central core taken from the NH3 emission peak. 
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Fig. 2. N2H + (l-0) spectra toward three positions of 
AFGL5142. The three positions labeled on each panel are, 
from top to bottom, CC (NH3 peak or central core), EC (peak of 
the eastern core), and WC (peak of the western core, see Fig. |4j». 



resolution of 0.42 km s to observe the N2H + (1-0) emission. 



ing 3C84. Data were cali brated and image d using the standard 
procedures in MIRIAD (ISault et al.l Il995h . We combined the 
data from both days of observations. The rms in the naturally 
weighted maps is ~ 0.7 mJy beam -1 in the continuum, and 
~ 55 mJy beam -1 per 0.42 km s -1 channel in the line data. 
The synthesized beam of the continuum image is 2'.'53 X l'/71, 
P.A.= 68.8°. For the line emission, which is more extended than 
the continuum, we applied a Mt)-taper function of 52 kA (3'.'5 in 
the image plane) to improve the signal-to-noise ratio and to re- 
cover the extended emission. The resulting synthesized beam is 
4'.'04 x 3'.'53, with P.A.= 79.1°. The continuum and N 2 H + emis- 
sions were cleaned using a box around the emitting region. 

3. Results 

3. 1 . Continuum Emission 

Figure Q] presents the continuum map at 3.2 mm. We detected 
a single millimeter peak associated with IRAS 05274+3345 
East (we could not resolve the two main 1.2 mm peaks 
MM-1 a nd MM-2, separate d by about 1", detected with the 
SMA by Zha ng et alj|2007b . The millimeter continuum emis- 
sion has a single compact peak component with a flux density 
of 107.1 + 1.6 mJy, surrounded by a more extended and flat- 
tened emission elongated toward the west and south-east di- 
rections. A 2D Gaussian fit to the compact component yields 
a deconvolved size of 2'.'9 x 273, P.A.= 16.5° (5200 AU at 
the distance of the source). The peak position is a{J2000) = 
05 h 30 m 48 s .01; 5(72000) = +33°47'5471. The total flux density 
taking into account the extended emission is 126.4+3.1 mJy, 
in very g ood agreement with the flux density at 88 GHz re- 
ported by iHunter et al.l d 1999h . which is ~125 mJy. Assuming 
that the du st emission is optic ally thin, and using the dust opac- 
ity law of iHildebrandl ( 1983 ). with a dust emissivity index of 
f3 — 1 (see Zhang ^talT 20071 for a description of its deriva- 
tion), and for a dust tempera ture of 45 K (based on the esti- 
mation of IHunter et al.l fl999f) the total mass of dust and gas 
from thermal continuum emission is 23 Mq (after correcting 
for the expected c ontribution, ~4 %, from the ionized gas; see 
Hunt er et a This mass is consistent with the mass esti- 

mated from the 1 m m SMA observations (~ 21 Mq). Although 
IZhang et"al] d2007l) report a mass of 50 Mq, the difference arises 
from the value adopted for the dust emissivity index, =1.5 in- 
stead of 1 . We note that there is a factor of 4 in the uncertainty 
of the mass, mainly due to uncertainties in the dust opacity and 
the dust emissivity index. 



Phase calibration was performed with the quasars 0530+135 
and 0555+398, with typical rms in the phases of 12° and 8°, 
respectively. The absolute position accuracy was estimated to 
be around 071. Flux and bandpass calibration was set by us- 



3.2. Molecular line emission: N 2 H + 

In Fig. [2] we show the N2H + (l-0) spectra, not corrected for 
the primary beam response, at some selected positions of 
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Fig.3. Top panel: N2H + (1-0) channel maps for the hyperfine F\F = 01 — * 12 (isolated line) toward AFGL 5142. Contour levels 
are -2, 2, 4, 6, 8, 10, and 12 times the rms noise of the maps, 0.055 Jy beam" 1 . Bottom panel: N2H + (1-0) channel maps for the 
hyperfine F\F = 21 -> 11, 23 -> 12, 22 -> 11, and FiF = 11 -> 10, 12 ^ 12, and 10 11. Contour levels are -4, -2, 2, 4, 6, 
8, 10, 12, and 14 times the rms of the map, 0.055 Jy beam" 1 . In all panels the cross marks the position of the 3.2 mm continuum 
source. The synthesized beam is shown in the bottom left corner. Velocities are referred to the FiF = 01 — > 12 hyperfine, with the 
systemic velocity i)lsr= _ 3 km s" 1 . 



AFGL 5 142, and in Table Q] we present the line parameters ob- 
tained from the fit of the hyperfine spectra toward these po- 
sitions. Figure [3] presents the velocity channel maps for the 
N 2 H + (l-0) emission toward AFGL5142. The channel with 
maximum intensity, for the hyperfine component FiF = 01 — > 



12, was found at the velocity v= -3.5 km s 1 (the systemic ve- 
locity of the cloud is i)lsr= - 3 km s" 1 ). 

The zero-order moment map integrated for all the hyperfine 
transitions is presented in Fig. [4^. We also show (in Fig. I^i), 
for comparison, th e NH3 (1,1) integrated intensity map from 
IZhang et alJ ((2002) overlaid with the cluster of infrared sources 
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Fig. 4. a) N2H + zero-order moment integrated for all the hy- 
perfine components of the (1-0) transition. Contours start 
at 3 %, increasing in steps of 10 % of the peak inten- 
sity, 6.01 Jybeam 'km s _1 . b) NH3 (1,1) zero-order moment 
dZhang et al.l2002l) . Contours start at 10 %, increasing in steps of 
10 % of the peak intensity, 0.202 Jy bearrr'km s _1 . The synthe- 
sized beams, 4'.'04 X 3'.'53 for N 2 H + and 3'.'58 x 2'.'64 for NH 3 , are 
shown in the bottom right corner. In both panels the blue cross 
marks the position of the millimeter source reported in this work. 
Red crosses: type I (protostars ); green crosses: type II (class II 
objects) from Qiu et al. (2008). Yellow and green tilted crosses 
are from lChen et al.l d2005l) and represent Class I luminous pro- 
tostars with M < 5 Mq, and Herbig Ae/Be stars or T Tauri 
stars, respectively. Whit e crosses are the 5 m illimeter peaks de- 
tected with the SMA by Zh ang et al.l (|2007), and the star marks 
the position of the IRAS 05274+3345. The western, eastern, and 
central cores are also labeled as WC, EC, and CC, respectively. 
The thick strai ght black lines represent the direction of outflows 
A, B, and C dZhang et al.ll2007l) . 



reported in rece nt studies dHunter et al.lll995t IChen et alJl2005l 
lOiu et al]l2008h in order to have a complete view of the stellar 
content in AFGL5142. The overall structure of the integrated 
N 2 H + (1-0) emission consists of two main cores, located to the 
west and to the east of the dust condensation (hereafter west- 
ern and eastern cores), surrounded by a more extended struc- 
ture with an X-shape morphology, roughly following the same 
morphology as the NH3 emission. As can be seen in Fig. |U 
the western core seems to contain at least one infrared source 
embedded in the dense gas, located close to the position of the 
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Fig. 5. a) First-order moment map (velocity) for the hyperfine 
F\F = 01 -> 12 line of N 2 H + (1-0) toward AFGL5142 (color 
scale), b) Second-order moment map (velocity dispersion) for 
the hyperfine F X F = 01 -> 12 line of N 2 H + (1-0) (color scale). 
In both figures black contours are the same as in Fig|4](a), show- 
ing the N2H + (1-0) emission, with contours starting at 3 %, and 
increasing in steps of 10 % of the peak intensity. Color scales 
are in km s -1 . The synthesized beam is shown in the bottom 
right corner of the image. The black cross marks the position of 
the 3.2 mm source reported in this work. Bl ue and red lines rep - 
resent the direction of outflows A, B, and C (Zhang et al. 2007), 
and the black dashed line in the top panel indicate the position- 
velocity cuts (see Fig Note that the second-order moment 
gives the velocity dispersion, and must be multiplied by the fac- 
tor 2 V2 In 2 =i 2.35 to convert to full width to half maximum. 

NH3 peak, and which displays typical characteristics of proto- 
stars (Che n et al. 1 12005b . whereas the infrared source lying close 
to the peak position of the eastern core has no associated infrared 
exce ss and has been classified as a Herbig Ae/Be or T Tauri 
star dChen et al.ll2005l) . However, there is no clear evidence that 
these infrared sources are really associated with the dense gas, 
and additi onally, since the sensitivity of Spitzer observations is 
~ 1 Mq (lOiu et al. 2008), and that of near-infrared studies is 
-0.2-0.4 Mq dChen et al.ll2005l) . we cannot discard the possibil- 
ity of more infrared sources embedded in the dense gas. A 2D 
Gaussian fit to the western and eastern cores yields deconvolved 
sizes of ~ 9'.'3 x 776 (~ 0.08 x 0.07 pc), and ~ 1074 x 873 
(~ 0.09 x 0.07 pc), respectively. Interestingly, toward the dust 
condensation where the intermediate/high-mass stars are deeply 
embedded in the strongest and compact NH3 core (hereafter cen- 
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tral core), the N2H + emission decreases drastically. However, it 
is important to note that the central core is detected in N2H + in 
some velocity channels, from 7 km s -1 to 8.7 km s _1 , whose 
peak positions are slightly shifted to the east (~ 2") from the 
dust and NH3 peaks. 

We used the hyperfine component F\F — 01 — > 12 to study 
the kinematics of the N2H + emission seen from Fig.[3]since this 
hyperfine component is strong and not blended with the other 
hyperfine components. In Fig. [5^ we show the first-order mo- 
ment map (intensity weighted mean olsrX an d m Fig- Ep the 
second-order map (intensity weighted velocity dispersion) for 
this hyperfine. 

The first-order moment map shows that the eastern and west- 
ern cores appear at different velocities (see also the channel maps 
shown in Fig. [3]). Toward the eastern core there is a small veloc- 
ity gradient in the east-west direction, from ~ -3.5 km s -1 to 
~ -1.8 km s _1 . On the other hand, the western core shows a 
velocity gradient, from ~ -4.7 km s -1 to ~ -2.6 km s -1 with 
increasing velocities from the southwest to the northeast. In fact, 
as can be seen in Fig.[5p the emission from the western core has 
very broad lines, with a velocity dispersion of ~ 0.8-1 km s _1 , 
which corresponds for a Gaussian line profile, to a full width at 
half maximum (FWHM — 2 V21n2 x cr„, where <x„ is the veloc- 
ity dispersion) of ~ 1.8-2.3 km s~' (corrected for instrumental 
resolution) at the eastern side of the western core. This region 
of broad emission coincides with the passage of outflow B of 
IZhang et"aL1 (120071) . The values around ~ 2 km s _1 are signifi- 
cantly higher than the thermal line width ~0.2-0.3 km s -1 (esti- 
mated for a kinetic temperature of 25 K and 45 K, respectively), 
indicative of N2H + having a significant contribution from non- 
thermal processes, such as turbulence injected by the molecular 
outflows, and/or global systematic motions. Regarding the east- 
ern core, the typical value found for the velocity dispersion is 
in the range ~ 0.3-0.5 km s _1 , corresponding to line widths of 
~ 0.2-1.1 km s~' (corrected for instrumental resolution). In ad- 
dition, there is a clumpy cavity structure toward the northwest 
clearly visible in the -2.6 and -2.2 km s -1 channel maps of 
Fig. [3] which appear redshifted by ~ 1 km s~' with respect to 
the systemic velocity. 

In Fig. [6^ we present the position-velocity (PV) plot made 
across the east- west direction (P.A.= 90°), encompassing the 
N2FF peaks of the western and eastern cores, which clearly 
shows a ring-like structure with the western core having very 
broad line widths. The line broadening, of ~ 3.5 km s , is also 
seen in the PV-plot made across the western core at P.A.= 15° 
(Fig.[6p). Additionally, we performed a third cut at P.A.= 140° 
(Fig. [(J), picking up the eastern side of the cavity and the eastern 
core. From this plot, it seems that there are two structures, a ve- 
locity gradient associated with the eastern core, of ~ 1.5 km s -1 
in a region of 7", and a curve structure associated with the cav- 
ity. Note that the curve structure seen in the PV-plot of Fig. [6J; 
suggests that the cavity is expanding. Overall it seems that both 
cuts, at P.A.= 90° and P.A.= 140°, reveal expanding motions, 
suggesting that these motions could be produced by a combina- 
tion of stellar winds, radiation, and/or molecular outflows from 
the central cluster. Hence, the AFGL 5 142 cluster seems to be in 
process of disrupting the natal cloud. However, we will not fur- 
ther discuss the kinematics of the region since it is not the aim 
of this paper. 
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Fig. 6. N 2 H + p-v plot for the F\F = 01 -» 12 hyperfine at a) 
PA=90° (along the eastern and western cores). Channel maps 
have been convolved with a beam 6" x 2", with P. A. perpendic- 
ular to the direction of the cut. Contours start at 10 %, increasing 
in steps of 5 % of the peak emission, 0.408 Jy beam -1 . The cen- 
tral position corresponds to the 3.2 mm source. Positive offsets 
are toward the east, b) PA= 15°. Channel maps have been con- 
volved with a beam 2" x 2". Contours start at 10 %, increasing 
in steps of 10 % of the peak emission, 0.601 JybeairT 1 . The 
central position is taken at the peak position of the western core, 
which is Ax = -6.5", Ay = -3" offset from the phase center. 
Positive offsets are toward the northeast, c) PA= 140°. Contours 
start at 20 %, and increase in steps of 10 % of the peak emission, 
0.376 Jy beam -1 . Positive offsets are toward the southeast. The 
central position corresponds to Ax = 3.9", Ay = -0.59" off- 
set from the phase center. In all panels, the vertical dashed line 
indicates the systemic velocity, ulsr =- 3 km s _1 . 

4. Analysis 

4. 1 . Column density maps 

We studied the chemical environment of AFGL 5 142 by an- 
alyzing the column density of several molecular species. We 
used the N 2 H + (1-0) (this work), NH 3 (1,1) and NH 3 (2,2) (from 
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Fig. 7. a) N2H + column density map. b) NH3 column density 
map. c) H I3 CO + column density map overlaid with the 3.2 mm 
continuum emission (contours). In all panels scale units are 
cm -2 . The synthesized beam, ~ 4", is shown in the bottom right 
corner of each panel. The white cross marks the position of the 
3.2 mm peak. 



2002), and HCO + (1-0) and H 13 CO + (1-0) (from 



Zhang et al 

Hunter et al. 19991) data. In order to properly compare the emis- 
sion of all the molecules, we convolved the N2H + , NH3, HCO + , 
and H 13 CO + channel maps to obtain a final circular beam of 
~ 4". We computed the column density maps by extracting the 
spectra for positions in a grid of 1" x 1". Using CLASS we fit- 
ted the hyperfine structure of each spectrum for N2H + (1-0) and 
NH 3 (1,1), and a single Gaussian for the NH 3 (2,2), HCO + (1-0) 
and H 13 CO + (1-0). We fitted only those positions with an inten- 
sity greater than 5cr for N2HL (1-0) and NH3 (1,1) in order to 
ensure we are detecting all the hyperfine components, whereas 
for the other molecular species we fitted the spectra with an in- 
tensity greater than Act. In the following sections we show the 
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Fig. 8. Color scale: N(NH3)/N(N2H + ) abundance ratio map over- 
laid with the N2H + emission (black contour), with contours start- 
ing at 3 %, and increasing in steps of 20 % of the peak intensity. 
The synthesized beam, ~ 4", is shown in the bottom right corner. 
The white cross marks the position of the 3.2 mm peak. 

Table 2. Average column density for each dense core 





W(NH3) 


N(N 2 H + ) 


W(H 1J CO+) 


A'(NH 3 )/A'(N2H + ) 


Core 


(cm" 2 ) 


(cm" 2 ) 


(cm -2 ) 




central 


2.6 x lO 15 


2.5 x 1() 12 


4.5 X 10 12 


1000 


eastern 


1.4 x 10 15 


1.1 x 10 13 




130 


western 


1.8 x 10 15 


1.8 x 10 13 




100 



method used to compute the column density for each molecule. 
The main results are summarized in Table [2] 



4.1.1. N 2 H + 

We used the hyperfine structure fitt ing program in CLASS 
dForveille. Guilloteau. & Lucasl [19 89) adopting the hyperfine 
frequencies given in ICaselli et al.l dl995h to determine the i>lsr, 
the intrinsic line widths, total optical depths, and excitation tem- 
peratures (T ex ) in each position of the grid. The value of 7 ex was 
derived assuming a filling factor of 1 . As shown in Table Q] the 
excitation temperature, 7 ex , derived from the hyperfine fits in the 
eastern core, is in the range ~ 14—16 K, while 7 ex has lower val- 
ues, around ~ 10-12 K and ~8-9 K, for the western and central 
cores, respectively. 

The results obtained from the fits indicate that the N2I-L 
emission is essentially optically thin for most of the region 
(ttot —0.3-0.6). However the optical depth in the western core 
reaches higher values, around ttot —1.5—2. We calculated the 
N2H + column density, corrected for the primary b e am resp onse, 
following the expression given in ICaselli et alJ d2002bl) . and 
approximating the partition function to Q ml (kT ex /hB) ~ 
0.4473 r ex , where k is the Boltzmann constant, 7 ex the exci- 
tation temperature, h is the Planck constant, and B is the rota- 
tional constant of N2H + . Figure [7^ shows the resulting N2H + 
column density map. Clearly, we found important variation in 
the N2H + column density in AFGL5142. The highest value of 
the N2H + column density, 3 x 10 13 crrT 2 , is reached in the west- 
ern core. Toward the eastern core we found values of the N2H + 
column density around ~ 1 x 10 13 crrT 2 , while the N2H + column 
density has the lowest values, around ~l-5xl0 12 cirT 2 , toward 
the millimeter condensation, i. e., toward the central core (see 
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Table 0. The values found in the western and eastern cores are 
in good agreement with the N?H + column den s ities r eported by 
iPirogov et all d2003l 120071) . and iFontani et ail |2006) for dense 
molecular cloud cores with massive stars and star clusters, and 
for a sample of high-mass protostellar candidates, respectively. 
In addition, the values of the N2H + column density found in the 
western and eastern cores of AFGL 5 142 are also consistent with 
those reported in recent studies conducted toward m assive star- 
forming regions observed with in terferometers (e. g., Pala u et alj 
l2007t Heuther & Henningl l2009). but those obtained in the cen- 
tral core are clearly below. It is important to emphasize the dif- 
ference of 1 order of magnitude in the N2H + column density 
in the central core compared to the western and eastern cores. 
The uncertainty in the N2H + column density is of the order of 
25-50 %, and it has been estimated taking into account the un- 
certainty in the output parameters of CLASS and the uncertainty 
in the calibration. 

4.1.2. NH 3 

The optical depth of the main line derived from the NH3 emis- 
sion is in the range ttot -1.5—2.5 toward the central and 
compact NH3 core containing the intermediate/high-mass stars, 
while the emission in the extended structure is optically thin. 

From the results of the fits of the NH 3 (1,1) and NH 3 (2,2) 
spectra we derived the rotational temperature (r iot ) and com- 
puted the N H 3 column density m ap f ollowing the proce dures 
described in iHo & Townesl d 19831) and lHariu et alj d 19931) (see 
also the appendix of iBusquet et alJl2009l for a description of the 
method). It is worth noting that the rotational temperature map 
obtained from this analysi s is similar to the N H 3 (2,2)/NH 3 (1,1) 
intensity map presented in Zhang et al.l d2002l) . 

The NH 3 column density map, corrected for the primary 
beam attenuation, is shown in Fig. [7J3. The map shows small 
variations along the entire cloud. Contrary to N2FL, the NH 3 col- 
umn density reaches the maximum value, around 3 x 10 15 cm" 2 , 
toward the central core, at the position of the millimeter conden- 
sation, and where we found the minimum N2H + column density 
(see Table|2]). The values found for the NH 3 column density to- 
ward AFGL 5 142 are similar to those obtained by iPalau et alj 
d2007l) in the massive star forming region IRAS 20293+3952. 
The uncertainty in the NH 3 column density is of the order 
of ~15-25 %, e stimated from the un certainty in the rotational 
temperature (see Busquet et al. 2009), the output parameters of 
CLASS, and the uncertainty in the calibration. 

4.1.3. H 13 CO + 

To derive the H 13 CO + column density, we assumed that 
H 13 CO + (1-0) emission is optically thin while HCO + (1-0) is 
optically thick. The H 13 CO + (1-0) optical depth has been esti- 
mated from the ratio of the H 13 CO + and HCO + integrated in- 
tensity emission obtained from Gaussian fits and excluding the 
high- velocity components of HCO + (1-0) coming from molecu- 
lar outflow emission. 

Assuming that all levels are populated according to the 
same excitation temperature r ex , and approximating the parti- 
tion function to kT ex /hB, with B being the rotational constant of 
the molecule, the column density of a linear molecule obtained 
from the transition J —* J — 1 is: 



N 



3k 1 exp(gp|g) 



where fi is the electric dipole momentum of the molecule, vio the 
frequency of the transition / = 1 — > 0, and to is the opacity at 
the line center. 

For the particular case of H 1 3 CO + (1-0) molecular transition, 
H = 3.88 Debye= 3.88 x 10" 18 (cgs), v 10 = 86.754 GHz, the 
column density can be written, in practical units, as 



JV(H 13 CO + ) 



2.56 x 10' 



1 



1 



-4.i6/r cx 



K 



TO 



Av 



km/s 



(2) 



T A V, 



(1) 



In Fig. [7J: we present the H 13 CO + column density map to- 
ward the central core of AFGL 5142. An obvious feature is 
that the H I3 CO + column density map is correlated with the 
NH 3 column density map, reaching the maximum value, around 
5 x 10 12 cm" 2 toward the millimeter condensation, but is anticor- 
related with the N2I-L column density. The values obtained for 
the H 13 CO + column density are similar t o those found toward re - 
gions of high-mass star formation (e. g.. lZinchenko eta3|2009). 
Taking into account the uncertainty in the output parameters of 
CLASS and the uncertainty in the calibration, we estimated an 
uncertainty in the H I3 CO + column density around ~30-45 %. 

4.2. The A/H 3 /A/ 2 H + abundance ratio map 

Figure [8] shows the NH 3 /N2H + abundance ratio map toward 
AFGL 5 142. First of all, it is important to remark that the N 2 H + 
spectra toward some positions of the central core have low 
signal-to-noise, so only a few positions with high enough signal- 
to-noise ratio were fitted properly, and hence the values obtained 
toward the central core, specially those toward the south western 
side of the millimeter condensation, are overestimated due to the 
interpolation process. A clear feature seen from the NH 3 /N2H + 
ratio map is that the highest values, ranging from ~400 up to 
1000, are reached close to the position of the central core (i. e., 
harboring a cluster of millimeter sources). We explicitly checked 
that such a high values are real, not produced by the interpola- 
tion process. Toward the western core, which could contain a 
low-mass YSO, we found values of the NH 3 /N2H + ratio around 
50-100, whereas in the eastern core, maybe containing a more 
evolved object, the NH 3 /N2H + ratio is around 50-100 toward 
the center of the core and increases up to 200 at the edges. In 
Table [2] we give the average value of the NH 3 /N2H + ratio for 
these three cores. The uncertainty estimated for the NH 3 /N2H + 
abun dance ratio is a round ~30-60 %. As already pointed out 
bv lZhangetaD d2002l) . the VLA NH 3 data recov er 65 % of the 
flux observed by the Effelsberg 100 m telescope (Estalella et al. 
[19931) . and the missing flux arises mainly from smooth structures 
with scales larger than 50". Unfortunately, we can not estimate 
the amount of flux filtered out by CARMA in N2FL due to the 
lack of short spacing information. However, since the CARMA 
N2H + and VLA NH 3 images are sensitive to similar angular 
scales, the NH 3 /N2H + ratio is not seriously affected by differ- 
ences in the missing flux. 



5. Chemical modeling 

5.1. The UCL_CHEM model 

We employed UCL_CHEM, a time- and depth- dependen t 
chemical model similar to that used in IViti & Williams] d~1999); 
IViti et all d200ll I2004I) to study the behavior of the NH 3 /N 2 H + 
abundance ratio. For simplicity, in this study we performed 
time-dependent single-point calculations. The model consists 
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Table 3. Initial gas-phase elemental abundances relative to hy- 
drogen nuclei. 



Atom Abundance 



H 


1.0 




He 


0.075 







4.45x10" 


-4 


C 


1.79x10" 


-4 


N 


8.52x10' 


-5 


S 


1.43x10" 


-6 


Mg 


5.12x10- 


-6 



Table 4. Parameters of the UCL_CHEM model 



Parameter 


cc a 


WC and EC" 


Depletion (%) 


99 


99 


Gas Density (cm -3 ) 


10" 


10 ? 


At 


41 


10 


yd 

max 


70 


25 



Notes. a CC: central core. b WC and EC: western and eastern cores. 
c At the end of phase I. d Maximum value reached during phase II. 



of a two-phase calculation: in phase I we follow the chemi- 
cal and dynamical evolution of a collapsing core. This phase 
starts from a fairly diffuse medium (n a 400 cm" 3 ) in neutral 
atomic form (apart from a fraction of hydrogen in H2) that un- 
d ergoes a collapse follow ing the free-fall collapse law described 
in Rawlin gs et al J (11992) until densities typical of massive dense 
cores (n ^ 10 5 - 10 7 cm" 3 ) are reached. During this time, atoms 
and molecules from the gas freeze o n to the dust gra i ns an d 
they hydrogenate where possible, as in IViti & Williams! <[[999). 
The degree of freeze out or depletion, defined in terms of the 
depletion of the CO molecule, is a free parameter (although 
linked with density, see Rawlings et al. 1992), which we explore 
in this study. During phase I we adopted a constant tempera- 
ture of 12 K, and we assumed a standard value for the cosmic 
ionization rate, £ = 3 x 10" 17 s _1 . The initial gas-phase ele- 
mental abundances relative to hydrogen nuclei adop ted in this 
work, based on the findings of Sofia & Mever (2001), are listed 
in Table [3] We included 127 gas-phase species and 42 surface 
species interacting in 18 69 chemical reaction s adapted from the 
UMIST 2006 database (IWoodall et alJl2007l) . During the col- 
lapse phase (phase I) we let the chemistry develop for 3 Myrs, 
regardless of when the final density is reached. 

In phase II we explored the chemical evolution of the core 
harboring a newly born massive star by simulating the presence 
of an infrared source at the center of the core, and by subject- 
ing the core to a time-dependent increase of the gas and dust 
temperature (assumed to be the same because of the high den- 
sities considered), which was performed differently depending 
on the core we were modeling. For the central core we adopted 
a maximum temperature of ~ 70 K, whereas the temperature 
toward the western and east ern cores is aro und ~ 25 K (both 
derived from NH3 data, see lZhang et alJ[2002h . During phase II 
we performed the same treatment f or the increasing temperature 
and evaporation from grains as in IViti et al] d2004l) . which in- 
cludes the experimental results on desorption from grains and 
the evaporation of icy mantles formed in star-forming regions 
dCollings et al.ll2003alfbl 120041) . In this case, the evaporation of 
a fraction of mantle species occur when t he temperature for a 
particular desorption event is reached (see lViti et al.l l2004 for a 
complete description). In phase II we followed the evolution of 
each core for 3 Myr. 



5.2. Results 

We computed a grid of models to investigate the differences ob- 
served in the N2H + column density and the NH3/N2H + abun- 
dance ratio. For the central core we followed the chemical evo- 
lution of a core of 0.02 pc in diameter, whereas for the eastern 
and western cores we adopted a size of 0.04 pc. We explored 
a range of densities, n 10 5 — 10 7 cm" 3 , and for a particular 
density we used different percentages of freeze-out at the end 
of phase I. From this analysis we found that the observed val- 
ues of the NH3/N2H + ratio can be reproduced by our model us- 
ing a high degree of depletion at the end of the collapse phase 
(phase I), and by adopting a density « a 10 6 cm" 3 for the cen- 
tral core and n a 10 5 cm" 3 for the western/eastern cores. It is 
worth noting that the model for the central core is supported by 
our estimate of the density from the dust continuum emission. 
Adopting the mass derived in this work (see Sect. 3.1) and a 
size of 6", we obtained a density of ~ 5 x 10 6 cm" 3 , consistent 
with the density adopted in the chemical model. In Table [4] we 
list the parameter choices for the two selected models (CC and 
WC/EC models), which can reproduce both the observed NH3 
and N2H + column densities and the NH3/N2H + abundance ratio. 
In the following sections, we discuss the main results obtained 
by our chemical model for phase I and phase II. 

5.2.1. Fractional abundances and the NH 3 /N 2 H + abundance 
ratio during phase I 

In Fig.|9](left panels) we present the fractional abundances (with 
respect to hydrogen nuclei) of some molecular species of inter- 
est (i. e., CO, N 2 , H+, NH 3 , N 2 H + , and HCO + ) together with the 
NH3/N2H + abundance ratio as a function of time for the central 
core (CC model) and the western/eastern cores (WC/EC model) 
during phase I. First of all, we note that the fractional abundances 
obtained for the central core (top panel) and the western/eastern 
cores (middle panel) do not show significant differences because 
they only differ in the core size and the final density reached at 
the end of this phase. As can be seen in Fig. [9] the fractional 
abundance of CO, N2, NH 3 , N 2 H + , and HCO + increases with 
time (with a more pronounced increase for the case of CO, N2, 
and NH3), reaching its maximum value at times t a 2.5 x 10 6 yr. 
After this time, the fractional abundances start to decrease be- 
cause molecules freeze out onto dust grains due to the low tem- 
perature (T - 12 K) and high density (n a 10 5 - 10 6 cm" 3 ) 
achieved in these cores. On the other hand, the ion keeps 
more or less a constant abundance until t a 2.5 x 10 6 yr, and 
then its abundance rises moderately. 

The primary formation route for N2H + is through the reac- 
tion between H and N2, which also produces H2. N2Ff + can then 
react with many atoms and molecules such as C, O, H2O or CO 
to form N2. Hence it is therefore clear that the abundance in the 
gas phase of these removal agents will be critical in determining 
the abundance of N2H + present. At high densitites, but before 
the effects of freeze-out dominate, all the molecular species that 
are important in the removal of N2H + and NH3 will be high. In 
particular, for both models the main formation route for N2H + 
is initially H2 + Nj — > N2H + + H, and then the chemical reac- 
tion Hi + N2— > N2H 4 - + H becomes important, while NH3 is 
mainly formed via dissociative recombination of NFT| . The lat- 
ter is formed via consecutive hydrogenation of N + (via reactions 
with molecular hydrogen). 

During phase I the NH3/N2FP abundance ratio behaves in 
the same way for both the central core and the western/eastern 
cores models, and only during the latest time steps it depends 
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Fig. 9. Left: Fractional abundances with respect to hydrogen as a function of time during phase I for the central core (top panel) 
and the western/eastern cores (middle panel). The NH3/N2H + abundance ratio for both the central core and the western/eastern core 
are also shown in the bottom panel. Right: Same as left panels but for phase II calculation. Note that the y-axis range according 
to the model. In the bottom right panel we also show the temperature as a function of time, indicated by the thin lines. The two 
shaded areas mark the range of values (400-1000 for the central core and 50-200 for the western/eastern cores) of the NH3/N2H + 
abundance ratio estimated from the observational data. 



on the adopted density (n - 10 5 crrT 3 for the western/eastern 
cores and n - 10 6 crrT 3 for the central core). In both cases, 
the NH3/N2H + ratio starts from low values and it reaches values 
around 1000 and even higher for / > 2.5 x 10 6 yr because CO de- 
pletion favors the formation of NH3 against N2H + in core centers 



at high densities. This suggests that long-lived starless cores are 
associated with high values of the NH3/N2H + abundance ratio. 
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time 


2V(NH 3 ) 


iV(N 2 H + ) 


NH 3 /N 2 H + 


Model 


(yr) 


(xlO 14 cm- 2 ) 


(xlO 11 cm" 2 ) 




CC 


1.9 x 10 4 


17 


38 


450 


cc 


4.5 x 10 s -5.3x10 s 


11-3.20 


0.80 


1400-400 


WC & EC 


1 x 10 4 -3xl0 6 


2.80-10 


56^11 


50-240 



5.2.2. Fractional abundances and the NH 3 /N 2 H + abundance 
ratio during phase II 

Figure [9] (right panels) shows the fractional abundances (with 
respect to hydrogen nuclei) as a function of time of CO, N2, 
HJ, NH3, N2H + , and HCO + derived from our chemical model 
during phase II for the central core (top panel) and the west- 
ern/eastern cores (middle panel). The fractional abundances in 
the central core clearly have a different behavior compared to 
that of the western and eastern cores as a consequence of the dif- 
ferent density, and especially the different temperature reached 
at the end of phase II. In the central core, the fractional abun- 
dance of CO, N2, N2H + , and HCO + is constant during the ini- 
tial stages, and once the temperature increases, their fractional 
abundances rise significantly. On the other hand, the fractional 
abundance of the molecular ion H3 decreases with time. This 
is of course due to the fact that is highly reactive at high 
densities. The abrupt increase of the fractional abundance takes 
place at the same time for all molecules, around f ^ 6.1 x 10 3 yr 
(corresponding to a temperature of T 21 K). This time cor- 
responds to the time when a fraction of weakly bound species 
evaporates from the grain mantle (see Table 1 in Viti et al. 
2004). The fractional abundance of N2H + reaches its maximum 
value 4.6 x 10" 11 (corresponding to a N2H + column density of 
~ 1.8 x 10 12 cm -2 ) at t =s 1.7 x 10 4 yr. The rapid destruction 
of N 2 H + (by the chemical reaction N 2 H + + CO -> HCO + + N 2 ) 
occurs at t ^ 2.3 x 10 4 yr, with a variation of the fractional abun- 
dance of one oder of magnitude, being the N2H + column density 
around ~ 8 x 10 1 1 cirr 2 , consistent with the value estimated from 
the observations. On the other hand, HCO + fractional abundance 
is constant at ~ 6 x 10 10 . Assuming [HCO + /H 13 CO + ]=40, we 
estimated an H 13 CO + column density of ~ 1.3 x 10 12 crrr 2 , in 
good agreement with the estimation from the observational data 
(see Table [2}. This chemical behaviour could be explained in 
terms of CO desorption from grain mantles, which causes a sub- 
stantial destruction of N2H + favoring the formation of HCO + , 
being this reaction the main removal mechanis m of N2H + for 
'standard' CO abunda nces of [CO/H2]- 10 4 (lJ0rgensen et alj 
2004; Le e et al.ll2004l) . In fact, as can be seen in Fig. [9] (top 
panel) the CO fractional abundance is ~ 10 4 , high enough to 
lead a significant destruction of N2H 4 ". 

Contrary to the case of N2H + and HCO + , the increase of 
the NH3 fractional abundance is slow, which indicates that gas 
phase chemistry dominates over pure evaporation. This is simply 
a consequence of the fact that most of the ammonia is still locked 
in water ice; in fact, as shown by the experiments by Collings et 
al. (2004), ammonia is released back into the gas phase only if 
temperatures of ~ 100-120 K are reached (Viti et al. 2004). The 
maximum value, of 4.1xl0~ 8 (or N(NH 3 )^ 3.2 x 10 15 cm -2 ), 
is reached at t 1.3 x 10 5 yr, and then it drops until it 
reaches a constant abundance of lxlO -9 (or a column density 
of ~ 8.2 x 10 13 cnT 2 ). This behavior, i.e., the decrease in the 
NH3 fractional abundance, takes place through reactions of the 
NH3 molecule with the ions C + and HCO + . In this situation, the 



C + fractional abundance increases at late stages, while the frac- 
tional abundance of HCO + is more or less constant, which is 
consistent with the fact that HCO + is produced from the destruc- 
tion of N2H + (see above) and destroyed through reactions with 
NH 3 . 

Regarding the western and eastern cores, which were mod- 
eled assuming the same core size (~ 0.04 pc) and same temper- 
ature (the maximum temperature is 25 K), we found that the 
fractional abundance of NH3 increases moderately with time, 
while the fractional abundances of CO, N 2 , N 2 H + , and HCO + 
rise considerably during the initial stages due to desorption ef- 
fects, similar to the case of the the central core model. At 
t ^ 4 x 10 3 yr, N2H + and HCO + have constant abundance, 
around ~ 2 x lO" 10 (N(N 2 H + )^ 4 x 10 12 cnr 2 ) and ~ 1.8 x lO" 9 
(N(HCO + )^ 3.6 x 10 13 cm 2 or N(H 13 CO + )=9xl0 n cm 2 , 
assuming [HCO + /H 13 CO + ]=40) for N 2 H + and HCO + , respec- 
tively, while the NH3 fractional abundance increases, up to ~ 
5.2 x 10- 8 (N(NH 3 )^ 1 x 10 15 cnr 2 ) at t =* 1.3 x 10 5 yr, and 
then it remains roughly constant too. For these cores, the values 
obtained from the chemical modeling are in agreement (within 
a factor of 2 in the case of N2H + , and a factor of 4 in the case 
of NH3) with the column densities reported from the observa- 
tional data. Finally, the CO fractional abundance in the west- 
ern/eastern cores is ~ 10 5 , significantly lower than the CO frac- 
tional abundance of the central core, indicating that the relatively 
low CO abundance in the western/eastern cores does not lead to 
a substantial destruction of N2H + . Therefore, the fraction of CO 
that will evaporate from grain mantles plays an important role 
in determining the fractional abundance of N2H + , and hence the 
NH3/N 2 H + abundance ratio. 

In Table [5] we show the values of the NH3 and N2H + col- 
umn densities, together with the NH3/N2H + abundance ratio for 
the two models that are in agreement with the observed values. 
In Fig. [9] (bottom right panel) we present the NH3/N2H + abun- 
dance ratio as a function of time for the central core and the 
western/eastern core obtained during phase II. The NH3/N2H + 
abundance ratio observed toward the central core, around ~400- 
1000, can be reproduced by our chemical model for times t 
10 4 yr and t =(4.5-5. 3)xl0 5 yr. For the time range t 10 4 - 
4.5 x 10 5 yr the model produces a higher NH3/N2H + abundance 
ratio, ~ 4000, slightly above the observed values for the cen- 
tral core. However we adopted the longer age as it is more re- 
alistic. For the western/eastern cores, the NH3/N2H + abundance 
ratio initially shows high values, then it decreases due to desorp- 
tion effects and finally reaches a constant value ~ 200 for typical 
ages of low-mass YSOs, at around t =s 10 5 - 10 6 yr. In addition, 
in Fig. [9] (bottom panel) we also show the temperature as a func- 
tion of time. For the central core, when the temperature is low, 
the NH 3 /N 2 H + ratio is high (> 10 4 ). Around T 21 K there 
is a clear drop of the NH3/N2H + ratio due to desorption effects. 
As temperature increases, the NH3/N2H + rises until it reaches a 
constant value, ~ 10 3 , at around T - 45 K. For the western and 
eastern cores the temperature varies from 12 to 25 K, producing 
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small variations on the NH3/N2H + ratio, and only for times in 
the range 10 3 -10 4 yr the NH3/N2H + ratio changes significantly, 
similar to the case of the central core. It is worth noting that in 
the model for the central core the visual extinction is Ay — 40, 
but typically the visual extinction in hot cores (i. e., embedded 
in the central core) is around ~ 100 or even higher. In order 
to evaluate the error in the NH3/N2H + abundance introduced by 
this difference we performed an additional model with a higher 
density, which gives a visual extinction of A v 100. In this 
situation, the NH3/N2H + ratio is affected by a difference of ~ 3- 
25 %.We also note that the abrupt drop of the NH 3 /N 2 H + ratio 
produced at early ages (f 5 x 10 3 yr) in both the central core 
and western/eastern cores is due to desorption effects of some 
molecules. 

Finally, we explored the fractional abundances of other 
molecular species and compared them with observations when 
available or made some predictions for future observations. We 
found that the fractional abundance of CH3CN obtained from 
the chemical model is in the range ~ 8 x 10~ 8 - 10~ 7 at 
t = (4.5 - 5.5) x 10 s y r , in agreement with the value reported 
by IZhang etal] d2007b . which is in the range (1 - 4) x 10~ 8 . 
Additionally, our model predicts a CO abundance of 10~ 4 in- 
dicating that CO is desorbed from grain mantles due to the 
relatively high temperature (around ~ 70 K) and the po wer- 
ful molecular outflows detected in CO (IZhang et al.l [2007). We 
found that molecules like HC3N, H 2 CO, C 2 H, and CS should be 
detectable toward the central core, since their fractional abun- 
dances are even higher than the NH3 ones. Concerning the west- 
ern and eastern cores, the CO fractional abundance, of 10~ 5 , is 
one order of magnitude lower than that of the central core. For 
these cores, molecules such as CS, C2H, and CH3CN are quite 
abundant and would be detectable but we would not expect to 
detect any H 2 CO and HC3N. 

In summary, the models that can reproduce the observed val- 
ues of the NH3/N2H + abundance ratio toward the central core 
as well as in the western and eastern cores, were obtained with 
a high depletion at the end of phase I. Hence, the differences 
between the central and the western/eastern cores seen in the 
NH3/N2H + abundance ratio seem to be mainly due to density 
and temperature effects. In particular, as shown by the chemi- 
cal modeling of the region, the temperature is a key parameter 
in determining the NH3/N2H + abundance ratio, and hence the 
abundance of each molecule. As temperature rises the fraction of 
species that will evaporate increase, and consequently the chem- 
istry evolves differently. 

6. Discussion 

The results obtained from the N2H + observation s, together with 
NH3 data from the literature (Zhang et al. 2002) show that there 
are three dense cores associated with the high-mass star-forming 
region AFGL5142. The central core, strong in NH3 but al- 
most devoid of N2H + , harbors a cluster of massive stars in their 
making. On the other hand, the western and eastern cores are 
stronger in N2H + than in NH3 emission. In the following sec- 
tion we discuss the possible reasons for the observed varia- 
tion of the NH3/N2H + abundance ratio among the cores found 
in AFGL5142. 

The central dense core in AFGL5142 is actively forming a 
cluster of massive stars, still in the accretion phase, which dis- 
plays the typical signposts of star formation (molecular outflow 
and maser emission). Concerning the western and eastern cores, 
the relatively high temperature, around 20 K, could be due to the 
presence of an embedded YSO and/or heating produced by the 



three molecular outflows associated with the central core, which 
are strongly affecting the surrounding dense gas (Zhan g~et al.l 
2002, |2007|) . Thus, the association of the western and eastern 
cores with an infrared source is not obvious. 

The NH3/N2H + abundance ratio map presented in Fig. [8] 
shows significant variations, with a ratio around ~50-100 toward 
the western and eastern cores, and a ratio, up to 1000, toward the 
central core, due to a significant drop in the N2H + abundance. 
The values found for the NH3/N2H + abundance ratio in the west- 
ern and eastern cores are similar to those found in previous 
studies for cores harb oring YSOs, for both low- and high-mass 
star-forming regions (Case lH et all l2002at iHotzel et alj l2004t 
Palau et al. 2007; Friesen et alj |2010h . These studies are all con- 
sistent with the fact that a high ratio (<300) seems to be associ- 
ated with starless cores and a low ratio, around 60-90, is found 
toward the YSOs, suggesting an anticorrelation between the the 
NH3/N2FU abundance ratio and the evolutionary stage. Thus, the 
high NH3/N2H + ratio toward the central core does not follow the 
anticorrelation between the NH3/N2H + ratio and the evolution- 
ary stage of the core like in the studies mentioned above. The 
noticeable increase of the NH3/N2H + ratio in the central core 
suggests a strong differentiation of the N2H + abundance between 
the central core and the western/eastern cores. Below we briefly 
investigate the origin of such a differentiation. 

UV radiation effects: 

As suggested by lOiu et al ] d2008l) . UV photons from 
IRAS 05274+334 are not likely to affect the dense gas, 
and the effects of the UV field from the embedded protostar(s) 
do not seem to be a major issue in determining the NH3 and 
N2FU abundance of the central core, since the high visual 
extinction, Ay — 40 mag or higher, prevents UV photons to 
penetrate the central core. 

Excitation effects: 

Even though we took into account the effects of opacity and 
different T ex in the calculation of the column densities, we 
consider whether the different excitation conditions among the 
dense cores may be the cause for the lower abundance of N2H + , 
i. e., the high NFWNp.H 4 " rati o. The critical density of N 2 H + is 
~ 10 5 crrT 3 (J0rgensen 2004), while the critical density of NH3 
is ~ 10 4 cm -3 (lHolll977l) . From the dust continuum emission 
we estimated a density of ~ 10 6 cirT 3 in the central core. In 
addition, the density obtained from single-dish observations at 
850 fim dJenness et al.l 1 19951) . adopting a core size of ~ 26", 
is ~ 4 x 10 5 cirT 3 . For these densities the transitions of both 
NH3 and N2FU are thermalized, and thus excitation effects do 
not seem appropriate to explain the N2H + emission drop in the 
central core. 

On the other hand, the hyperfine structure method in CLASS 
assumes equal excitation tem perature for all the hyperfine com- 
ponents. Danie feTaT] (|2006) show that this assumption does not 
hold for the case of high opacities due to radiative processes. 
Since in AFGL 5 142 we derived an opacity of ttot —0.3-0.6 for 
the N2FU molecule, the problem reported bv lDaniel et alJ {2006) 
is not likely affecting our result. 

Physical and chemical effects: 

The high values, up to 1000, of the NH3/N2FU abundance ra- 
tio found in the central core of AFGL 5 142 can be reproduced 
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by our chemical model for high densities (n =2 10 6 cm -3 ) 
and high temperaures (T ^ 70 K). The chemical modeling 
performed in Sect. 5 indicates that both density and tempera- 
ture play an important role in determining the molecular abun- 
dance of NH3 and N2H + , and hence their ratio. The cen- 
tral core has a higher temperature and density than the west- 
ern and eastern cores. Thus, a different chemistry can develop 
due to CO evaporation from the grain mantles. The CO des- 
orption in the central core of AFGL5142 leads to the destruc- 
tion of N2H + , and hence the NH3/N2H + abundance ratio in- 
creases considerably relative to the value found in the west- 
ern and eastern cores. This is supported by the fact that 13 CO 
is not frozen out in the central core, w hereas it is faintly de- 
tected in the western and eastern cores (Zhan gltaLll2007l) . As 
pointed out in Sect 5.2.2, Collin gset al.l ( 20041) find experimen- 
tal evidences that the desorption of NH3 from grain mantles 
takes place at a temperature of ~ 120 K. Therefore, since in 
our model we assumed a maximum temperature of 70 K, the 
high NH3/N2H 4 " abundance ratio in the central core is mainly 
a consequence of the destruction of N2H + by CO rather than 
an enhancement of NH3 in the gas phase. The disappearance 
of N2H + from the gas phase has been reported by several au- 
thors, both in low- and high-mass star-forming regions. In the 
low-mass regime, there are several reports of central N2H + de- 
pletion in Class 0/Class I protostars ( e.g., IRAM 0419 1 + 1522 : 
lBelloche& Andrei I2004t V LA 1623: iDi Francesco et all [2004; 
L483 mm: lj0rgensenll2OO4[: Barna rd lc: iMatthews et alJ l2006; 
BHR 7 1 IRS1: IChen et all l2008h . Regarding the high-mass 
regime, Pirog ov et al. ( 20031 12007) observed with a single-dish 
telescope a sample of dense cores associated with massive stars 
and star clusters containing IRAS point sources, and find that, for 
most of the sources, there is a decrease of the N2H + abundance 
toward the dust column density peak. Furthermore, interferomet- 
ric observations of N2H + toward the high-mass s tar-forming re- 
gions IRAS 23033+5951 dReid & Matthewsll2008l) also reveal a 
significant destruction of N2H + . Therefore, the high ratio mea- 
sured toward the central core of AFGL5142 seems to be due to 
the rapid destruction of N2H + when CO is released from grain 
mantles as seems to be the case for the other regions cited above. 

Evolutionary effects: 

Although in our chemical model we did not include molecular 
outflows, they may strongly affect the surrounding dense gas, 
an d hence the NH3/N 2H + abundance ratio. As already discussed 
by IChen et a"D d2008l) . who propose three stages of the interac- 
tion between the N2H + emission and molecular outflows, based 
on the morphology of the N2H + dense gas and the jet/outflow 
emission. Here we suggests a qualitative picture of the evolu- 
tion of the NH3/N2H + abundance ratio. First, in starless cores 
with high densities most molecules are highly depleted, favoring 
the formation of NH3 and N2H + . It seems that the CO deple- 
tion favors the formation of NH3 against N2H + in core centers 
at high densities. Thus, the NH3/N2H + ratio reaches high val- 
ues, consistent with t he results found in low-mass starless co res 
(ICaselli et al.ll2002at iHotzel et alJ[2004t iFriesen et alJl2010h as 
well as for starless cores nearby massive stars, like in the 
high- mass star-forming regions IRAS 20293+3952 dPalau et alJ 
2007). The high NH3/N2H 4 " ratio in starless cores has been re- 
produced by our chemical model shown in Fig|9] (left panels). 
After this phase, once the star has formed, the N H3/N2H 4 " ra- 
tio decreases, as found in the low-mass regim e dCaselli et alJ 
I2002at IHotzel et alJ l2004t IFriesen et alJ 1201 Ol) . an d in YSOs 
close to the massive star (e.g., IRAS 20293+3952: IPaTau et alJ 



12001 IRAS 001 17+612: Busquet et al. in prep.). In this situa- 
tion, CO is less depleted, and N2H + molecules in the envelope 
are entrained by the molecular outflow, being the N2H + emis- 
sion elongated in the d irection of the mo lecular outflow, as found 
in IRAS 20293+3 952 dPalau et alj2007l) and IRAS 001 17+6412 
dPalau et alj|2010l) . Finally, the next evolutionary stage involves 
the release of CO by powerful molecular outflows, which leads 
to the destruction of N2H + and forms large holes in the envelope. 
Then, the NH3/N2H + abundance ratio raises considerably, like 
in the central core of AFGL5142. 

Therefore, either molecular outflows and den- 
sity/temperature effects seem to be key pieces in the deter- 
mination of the NH3/N2H + abundance ratio. In order to quantify 
the importance of the temperature in the NH3/N2H + ratio, we 
compared the results obtained in AFG1 5142 with the low-mass 
star-forming region L483. This region seems to be in the same 
stage of disruption of the envelope as AFGL5142, as the N2H+ 
shows two clumps on both sides of the YSO (see J0r gensenl 
12004 iFuller & Woottenfe OOO). being then in the later evolution- 
ary stage of the scenario proposed above. However, the NH3 
emission in L483 is also showing a hole close at the position 
of the YSO, suggesting that possibly the NH3/N2H + abundance 
ratio is much lower than the values found in AFGL 5 142, around 
400-1000. This seems to indicate that the effects of molecular 
outflows alone cannot explain the large values derived in 
AFGL 5 142, although to draw a firm conclusion this should be 
modeled. Thus, the high temperature reached in the central core 
due to the presence of hot cores seems to affect significantly the 
NH 3 /N 2 H + abundance ratio, indicating that the NH3/N 2 H + ratio 
behaves differently than in low-mass star-forming regions. 

7. Conclusions 

We used CARMA to observe the 3.2 mm continuum emis- 
sion and the N2H + molecular line emission toward the mas- 
sive protostellar cluster AFGL 5 142. The N2H + dense gas emis- 
sion, which mimics the morphology of NH3, shows two main 
cores, one to the west (the western core) of the dust condensa- 
tion of ~ 0.08 pc of size, and another core to the east (the east- 
ern core), of ~ 0.09 pc. Interestingly, toward the dust conden- 
sation, of 23 Mq, the N2H + emission drops significantly, 
whereas the compact NH3 core (i. e., the central core) peaks at 
this position. The analysis of the N2H + column density and the 
NH3/N2H + abundance ratio indicates a strong differentiation in 
the N2H + abundance. While we found low values, around ~50- 
100, of the NH3/N2H + abundance ratio associated with the west- 
ern and eastern cores, the NH5/N2H + ratio rises significantly, 
up to 1000, toward the central core, in which the formation of 
several massive stars is taking place. We performed a chemi- 
cal modeling of the region using the time-dependent chemical 
model UCL_CHEM. The chemical model shows that a high 
NH3/N2H + abundance ratio can be reproduced in three differ- 
ent situations: i) in long-lived starless cores, ii) in high density 
cores with low temperature, ~20-25 K, and containing YSOs, 
and Hi) in high density cores with YSOs and high tempera- 
ture. For the central core we found that such high value of the 
NH3/N2H 4 " abundance ratio could be explained with a density of 
~ 10 6 cm -3 and a maximum temperature of 70 K, producing a 
substantial destruction of the N2H + molecule by CO, which is 
desorbed from the grains mantles due to the high temperature 
and possibly the molecular outflows associated with the massive 
YSOs embedded in the central core. On the other hand, to repro- 
duce the observed NH3/N2H + abundance ratio toward the west- 
ern and eastern cores, we used a lower density, of ~ 10 5 cm -3 , 
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and temperature (maximum temperature reached is 25 K). These 
results suggests that, in addition to the core evolution, the phys- 
ical properties of the core, such as density and temperature, play 
an important role in determining the NH3/N2H + abundance ra- 
tio. In conclusion, in AFG15142 we found that the NH 3 /N 2 H + 
abundance ratio behaves differently to the low-mass case mainly 
due to the high temperature reached in hot cores. 
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